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Abstract. In this survey paper, I first review the history of Bernouiii 
numbers, then examine the modern definition of Bernouiii numbers and 
the appearance of Bernoulh numbers in expansion of functions. I revisit 
some properties of Bernoulli numbers and the history of the computation 
of big Bernoulli numbers. 



1. Introduction 

Two thousand years ago, Greek mathematician Pythagoras first noted 
about triangle numbers which is l + 2 + 3 + -- -+n. Archimedes found out 

(1) 1^ + 2^ + 32 + • • • + = in(n + l)(2n + 1). 

Later in the fifth century, Indian mathematician Aryabhata proposed 

(2) 1^ + 2^ + 33 + • • • + = [^n(n + 1)]^ 

which Jacobi gave the first vigorous proof in 1834. It is not until five hundred 
years later that Arabian mathematician Al-Khwarizm showed 

(3) 1"^ + 2"^ + 3"^ + • • • + = — n(n + l)(2n + + 3n - 1). 

30 

n-l 

Studies the more generalized formula for ^ fc"" for any natural number r was 

only carried out in the last few centuries. Among them, the investigation of 
Bernoulli numbers is much significant. 

In this paper, I present an elementary examination of the development of 
Bernoulli numbers and a concise review of its appearance in expansions of 
various functions . I also aim to explore its properties and its applications 
in other fields of mathematics 



2. Bernoulli Numbers 

Swiss mathematician Jakob Bernoulli (1654-1705)once claimed that in- 
stead of laboring for hours to get a sum of powers, he only used several 
minutes to calculate sum of powers such as 1^'^ + 2^'^ + 3^'^ + • • • + 1000^*^ = 
91, 409, 924, 241, 424, 243, 424, 241, 924, 242, 500. Obviously, he had used a 
summation formula, knowing the first 10 Bernoulli numbers [2j. 
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It was already known to Jakob Bernoulli that 



n— 1 



(4) = ^"(^ - 1) = - 



k=l 

n-1 



1 

2- 2' 
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(5) ^A;2 = -n(n-l)(2n-l) = 3^''^ " + 

k=l 

(6) x:^'=z-v-i)' 



jfc=i 



1 ^ 1 o 1 2 

_n - -n + -n 



n-l 



V A:"^ = ^n(n - l)(2ra - l){3n^ - 3n - 1) = ^n^ - Jn^ + 

oU ^ 



k=l 



1 3 1 

(7) + ttW - —n 

^ ^ 3 30 

1 6 1 5 

= TT?^ - 7:n + 



n—l ^ 

^k^ = Y^n\2n'^ - 2n - l)(n - 1)^ 



k=l 

(8) +^»^-^»^ 



n-l 

More generally, the can be written in the form of 

k=l 



n-l r „ , 

= V ^' „r-fc+l 

^ ^ A;! (r- A; + l)! 

ifc=l fc=0 ^ ^ 

,0, ^0 ^ 5l . ^ i?2 ^ 



where the numbers which are independent of r and called Bernoulli's 
numbers. The first few Bernoulli numbers Bn are 
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-Bio 
Bl2 
Bi4 

Bie 
Bis 
B20 
B22 
B2A 
B26 
B28 
B30 
B32 

^34 
B36 

B38 
B40 



_5_ 
66 

691 



2,730 

7 
6 

3,617 



510 

43, 867 

798 

174,611 

330 
854,513 

138 

236,364,091 
2,730 
8,553,103 
6 

23,749,461,029 
870 

8,615,841,276,005 
14,322 
7, 709,321,041,217 
510 

2,577, 687,858,367 
6 

26, 315, 271, 553, 053, 477, 373 
1,919,190 
2,929,993,913,841,559 
6 

261, 082, 718, 496, 449, 122, 051 



13,530 
with 

(10) B2n+1 = 

for all positive integer n. More numbers are given in p]. 

3. Expansion of Usual Functions 

An equivalent definition of the Bernoulli's numbers is obtained from the 
series expansion of the identity 



n=0 



This leads to 



12) + - = -( + 1) = = -coth(-). 

' e^-1 2 2^e^-l ^ 2 - 1 2 ^2^ 
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It can be rewritten as 

(13) ^cothi-) = y^^. 

n=0 ^ ' 

If substitute x with 2ix, then it gives 

(14) -^ot{x) =±{-1)-?^^ 

for X G [— 7r,7r]. Thus the fohowing expansions are obtained, 

2B2„(2x)2"-i 



(15) coth{x) = ^ 



(16) cot{x) = ^(-1 



(2n)! 

n=0 ^ ' 

2BoJ2x?"-^ 



(2n)\ 

n=0 ^ ^ 

Now it's also possible to find the expansion for tanh{x) and tan{x). As I 
observe that 

cosh(2x) cosh(z) 



2coth(2x) — cothix) — ^ 

sinh{2x) sinh[z) 

cosh^ix) + sinh'^ix) cosh(x) 

~i:r~\ JT~\ ■ . / N = tanh{x) 

smn[x)cosn[x) sinn[x) 



I can have 



nr^ , u ^ ^ ^(4" - l)i?2n(2x)^"-i 
(17) tanh{x) = 



n=l 

and then, 



n=l ^ 

both for X € (— ^, ^) Bernoulh's numbers also occur in the expansions of 
other classical functions such as csc{x) ,csch{x) ,ln\sin{x)\,ln\cos{x)\,ln\tan{x)\, 
■ 1/ and etc. 

sinh(x) ' 

Another intriguing fact is that zeta function C(2fc) for any natural num- 
ber k can also be expressed in Bernoulli numbers (4j 

n=l 

However, there is no similar expression known for (^(2A; + 1). Rearrange the 
equation, I get 

(-l)^-i2(2fc)! ^ 1 (-l)^-i2(2fc)! ^,^,^, 

n=l 

This expression gives rise to an approximation of Bernoulli numbers. When 
k becomes large, 

(21) cm ~ 1, 
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while with Stirhng's formula 



(22) (2A;)! ~ {2kf''e~^''V4^, 
so I have 

(23) B2k ^ {-l)>'-H{—f''V^. 

ne 

4. Properties of Bernoulli Numbers 
The Bernoulli numbers are given by the double sum 

(^^> ^-E^iiD-D-C)'- 

k=0 r=0 ^ ^ 

where (^) is a binomial coefficient. They also satisfy the following interesting 
summations 0101 



n-l 



(25) EU"'^^ = °' 



fc=0 



n 



n 



(26) yA'lV]Bek = 2n + l, 



^\ 6k 

k=0 



I — n V / 



k=0 

The famous Clausen-von Staudt's theorem regarding Bernoulli numbers' 
fractional part was published by Karl von Staudt (1798-1867) and Thomas 
Clausen (1801-1885) independently in 1840. It allows to compute easily the 
fractional part of Bernoulli's numbers and thus also permits to compute the 
denominator of those numbers. It says, the value -62^1 added to the sum of 
the inverse of prime numbers p such that {p — 1) divides 2k, is an integer 
[7]. In other words, 

(28) -B2k= y -{modi) 



(p-l)|2fe^ 



For example, k 



-5i6 = 2^ - = - + - + - + — = — modi 
^-^ p 2 6 b n 510 

(p-l)|16^ 

463, 

=^ i3i6 = [modi). 

510^ ' 

One of the easy consequences of the Staudt's theorem is that for every prime 
numbers k of the form 3n + 1 

(29) B2k = \{modl). 

o 

This is because that p—1 divides 2k = 2(3n + 1) only if p — 1 is one of 1, 2, 
3n + 1, 6n + 2, that is p is one of 2, 3, 3n + 2, 6n + 3. But 6n + 3 is divisible 
by 3 and 3n + 2 is divisible by 2 because 3n+l is prime so the only primes 



6 



CONG LIN 



p candidates are 2 and 3. The first primes of the form 3n+l are 7, 13, 19, 
31, 37, 43, 61, 67, 73, 79, 97, 103... hence I have 

-Bl4 = B26 = i?38 = Bq2 = -B74 = ^86 = B122 

= ^134 = ^146 = ^158 = -B194 = B20Q = ■ ■ ■ = -(modi) 

6 

Staudt's theorem is very useful and significant in the sense that it permits 
to compute exactly a Bernoulli's number as soon as there is a sufficiently 
good approximation of it. 

Besides its applications in series expansions, Bernoulli numbers are also 
widely used in differential topology, mathematical analysis, and number 
theory. Interestingly, it is also related to the famous Fermat's last theorem. 

Fermat's Last Theorem states 

has no non-zero integer solutions for n > 2. Ever since Fermat expressed 
this result around 1630, generations of mathematicians have dived enthusi- 
astically into the pursuit of a vigorous proof . 

A breakthrough was made in 1850 by Ernst Kummer (1810-1893) when 
he proved Fermat's theorem for n=p, whenever p is a regular prime. Kum- 
mer gave the beautiful regularity criterion: 

p is a regular prime if and only if p does not divide the numerator of 
B2, B^, Sp_3. 

He showed that all primes before 37 where regular, hence Fermat's theo- 
rem was proved for those primes. 37 is the first non regular prime because 
it divides the numerator of 

7709321041217 37208360028141 

B^o = = 

510 510 

The next irregular primes (less than 300) are 

59, 67, 101, 103, 131, 149, 157, 233, 257, 263, 271, 283, 293, ... 
For example, 157 divides the numerators of Bq2 and -Biio- 

By applying arithmetical properties of Bernoulli's numbers, Johann Lud- 
wig Jensen (1859-1925) proved in 1915 that the number of irregular primes 
is infinite |Hj. 

5. Computation of Bernoulli Numbers 

Bernoulli himself calculated the numbers up to -Bio- Later, Euler worked 
up to B3Q. One century later, Adams made the computa of all Bernoulli's 
numbers up to -B124 0. In 1996, Simon Plouffe and Greg J. Fee computed 
-6200,000) and this huge number has about 800,000 digits. In July 10th 2002, 
they improved the record to -6750,000 which has 3,391,993 digits by a 21 
hours computation on their personal computer ^10,. The method is based 
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on the relation between zeta function and Bernoulli numbers, which allow 
a direct computation of the target number without the need of calculating 
the previous numbers. 
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